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We discuss differences between the exact S-matrix for scat- 
tering on serial structures and a known factorized expres- 
sion constructed of single-element S-matrices. As an illus- 
tration, we use an exactly solvable model of a quantum wire 
with two point impurities. 
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The core of the scattering-operator concept is that it relates asymp- 
totic states far away from the region where the interaction affects the 
motion substantially It happens often, however, that the interaction 
is localized to a certain domain Q ; since solutions to the free equation 
are usually known, we may rephrase then the scattering problem as a 
map between the solutions at the boundary of Q which can be con- 
tinued (in some direction) into the corresponding asymptotic states. If 
we want to keep a terminological distinction, it is more appropriate to 
speak about the prescattering operator in this case. 

Such a "finite-distance" scattering is particularly useful in situa- 
tions when the interaction support is a union of a finite number of 
domains fl, ; it can help to solve the full problem by means of scatter- 
ing on the "components". This subject has become actual in connec- 
tion with recent studies of quantum-wire superlattices [1-13] in which 
a number of elements, usually of the same type, is arranged into a 
serial structure. The physically relevant quantity is the conductivity 
which is related directly to the electron scattering in the superlattice 
by Landauer's formula. 

There are several ways how to derive the S-matrix components of 
the serial structure, i.e., transmission and reflection amplitudes from 
analogous quantities of a single element. The above sketched obser- 
vation applies directly if the number of linearly independent solutions 
entering and leaving each component scatterer is finite. This is the 
case, e.g., when the system in question is a graph || £§]; it may not be 
easy to find the solution interconnecting the scatterers if an external 
field is applied, but the algebraic part of the problem is well established 

However, quantum mechanics lures always around ready to show 
who is the master of our physical world. In general, diferent compo- 
nents of the wave function, say, different transverse modes in a quantum 
wire, do not cease to be correlated even out of the support of the in- 
teraction. Here the difference between the scattering and prescattering 
operator shows, because the sets of states the latter maps onto each 
other are larger; in addition to true asymptotic states which live eter- 
nally they contain also such that die out when the distance from the 
interaction region increases. 



2 





y 


d 


• 




• 


X 







Figure 1: A pair of scatteres in the strip. 



It would be certainly worth to formulate properly relations between 
this approach, which is inherently time-independent, and the rigorous 
scattering theory |R| |T7|. Our aim here, however, is more practical and 
connected with the mentioned studies of quantum-wire superlattices. 
With few exceptions their authors include the evanescent states into the 
iterative procedures of computing the S-matrix, so the result suffers 
no theoretical defect. The weak point of all numerical calculations 
stems from the necessity to restrict the used family of states to a finite 
number; the stability aspect is usually handled by a vague observation 
that the involved series converge fast enough. In this letter we want 
to discuss this problem in more detail; we are going to derive explicit 
error estimates in a solvable model in which the individual scatterers 
are point impurities in a straight strip. 



1 S— matrix factorization 

Elements of a one-dimensional superlattice are linearly arranged and 
most authors use various transfer-matrix modifications connecting so- 
lutions to the left and right of a given scatterer [1,2,5,9-11]; some admit 
numerical problems due to a fast growth of higher evanescent modes. 
It has been argued recently |L2] that the approach, which we called 



pre-S-matrix, offers a more stable scheme relating instead the ingoing 
and outgoing waves. 

Consider therefore a pair of scatterers Si, S 2 as sketched on Fig- 
ure 1, where A + is the family of solutions entering Si from the left 
etc.; in particular, if the corresponding external motion is free, A + is 
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Figure 2: A pair of point scatteres in the strip. 

a column vector (a^ e lkeX ^j with ke being the channel momenta. The 
pre-S-matrices can be written in terms of the reflection and transmis- 
sion parts, 

a) 

(1) 

where the tilded quantities corresponding to the passage from the right 
to the left are obtained from Rj, Tj by mirror transformation. The 
above form is convenient because the numbers of involved modes at the 
two sides of Sj may be different. Let us stress that the pre-S-matrices 
in this setting need not be unitary. The evanescent modes — if included 
- do not contribute to the probability current, and the propagating 
ones enter with different velocities: if the corresponding part of Sj 
should be unitary, we have to multiply the (n, m)-th amplitudes by 

( km I fori ) ^ • 

Our aim is to express D + from given A + , C_ . The relations 
(fj) yield a set of four operator equation which can be solved if we 
identify the left and right traveling solutions between the scatterers, 
D- = A- and C + = B + . Using the identity 

I + X(I-YX)~ 1 Y = (I-XY)- 1 

which is valid whenever the inverses make sense, we find by a straight- 
forward algebra that the combined pre-S-matrix is 

/ J R 1 +f 1 (/- J R 2J R 1 )- 1 i? 2 T 1 f^I-RiRj- 1 ^ 
Sn= | . (2) 

V T^I-R^y 1 ^ R 2 +T 2 (I-R 1 R 2 )- 1 R 1 f 2 
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The blocks of this operator acquire an illustrative meaning if we expand 
the inverses into geometric series ||. Sometimes S2 is obtained by 
shifting and/or mirroring another scatterer, then one has to find the 
operator which represents this transformation. This is equivalent to 
finding the general solution in the intermediate region between S\ and 
S2 which may not be easy if the motion is not free there ||15|| . 

The operator (0) describes the combined scatterer exactly if the 
space into which Si and S^ 1 map is large enough to accomodate all 
generalized eigenvectors of the corresponding Schrodinger equation. A 
truncation induces an error which in general is not easy to estimate; 
we shall do that below in an explicitly solvable model. 



2 Point— interaction scattering in a strip 

The model describes a nonrelativistic particle confined to straight strip 
of a width d with the hard-wall boundary subject to point perturba- 
tions simulating natural or artificial impurities. Such a system has nu- 
merous interesting features which will be described in detail elsewhere 
rS |; here we restrict ourselves to a basic information. 

For simplicity we set d := it and h 2 /2m = 1 , so that in the ab- 
sence of the impurities the motion is governed by the Hamiltonian 
H := —A ; the wavefunction is supposed to satisfy the Dirichlet 
boundary conditions ip(x,0) = ip(x,ir) = for any x. Point in- 
teractions situated at dj = (aj,bj) , j = 1, . . . , J can be introduced in 
the standard way |19|, Sec. 1. 5]: they are determined by the boundary 
conditions 

Li(ip, dj) + 271-ctjLoO, dj) = 0, j = 1, . . . , J, 
relating the generalized boundary values 

1p (x) r 

L (ip,d) := lim - — — — , Li(ip,d) := lim ip(x)— L (ifj, a) In \x—a\ 

\x-a\->0 m \X — a\ |x-a|-+0 L 

where otj are the (rescaled) coupling constants; the free Hamiltonian 
corresponds to a,- = 00 , j = 1, . . . , J . 
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/ 1 , \ i 00 

£(a: z) := lim [Go(a, x;z) In If— all = — 

SV ' ; |*-3|-o V V ' ' ; 2vr 1 7 7T ^ 



Solvability of the model stems from the fact that one can compute 
the corresponding resolvent kernel by means of the Krein's formula; all 
spectral information is contained then in the J x J matrix A(z) given 
by 

A(z) = 5 jm (otj - z)) - (1 -5 jm ) G (%, a m ; z) (3) 

where Go(xi,x 2 ]z) := ± £^° =1 g ~H^) ' 2l sin(nyi) sin(m/ 2 ) is the free 
Green's function, fc n (z) := \fz-ri 2 , and 

1 1 /^sin 2 (n6) 1 

fc n (z) 2in 
(4) 

We suppose everywhere that the energy stays away of the thresholds, 
\fz 7^ 1,2,.... The knowledge of the resolvent allows us, in particu- 
lar, to solve the scattering problem. The reflection and transmission 
amplitudes are expressed through the inverse of A(z) as 

(5) 

/ % . i , / s -iSm(mbj)sm(nb k ) _ i(k n ._ k ni) 

tnm(z) = 8 nm + - A (z)jk S , \ e • 

The mirrored quantities are obtained by changing each perturbation 
longitudinal coordinate a,j to — a,- ; together they satisfy the unitarity 
condition 



m=l m=l 



km(t nm t sm -\-T nm T srn ^ dns^n j ^ ' k m {t nm T sm -\- T nm t srr ij , 



where the summation runs over the open channels, [•] being the integer 
part. 



3 Comparison of the S-matrices 

From now on we shall consider a pair of point interactions, J = 2 , 
with ai — and a 2 = a . For the sake of brevity we denote jj := 
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aij—£(aj; z) and G := Go(«i, a 2 ; z) , so the matrix A(z) may be written 
as ( 2q 7^ ) • ^ or further purposes we introduce also the truncated 



resolvent 



G 1 = Gl(a u a 2 - z) := 1 f gM^^W ^ (6) 

among the M modes involved, A" := [\fz\ are propagating while 
the rest corresponds to evanescent channels. The one-perturbation 
reflection and transmission are 

u) _ i sin(Tt6j) sinjmbj) i(kn+km)a , {j) _ , i sin(r^) smjmbj) i(fc „- fcm)a ,. 

' nm 7 c 5 l mn u nm~T , e 

(7) 

for j = 1,2. Suppose that we employ an M component Ansatz for 
wavefunctions in the intermediate region, < x < a . Using the 
definition (|6|), we find 

f 1 - f«r^ = S nm - i G* Sin( " 6l) Sin(m&2) e*»" . 

This M x M matrix is explicitly invertible, 

v /nm vr 7i72-(G t )" ! k m 

The relations (0) and (H) allow us to express the composed transmis- 
sion amplitude as the lower left element of the matrix (|2|); after a 
straightforward computation using the definition (^j), we arrive at the 
formula 

i k~ l ( 

4m = 5 nm + ^TrHvI 1 72 sin(nfci) sin(m6i) + G* sin(n& 2 ) sin(m&i) e lfc " a 

vr 7i72-(<S*) i I 

(9) 

+G* sin(nfci) sin(m& 2 ) e~ ikma + 71 sin(n6 2 ) sin(m6 2 ) e i(fcn_fcm)a 

which differs from the exact expression (|5|) just by replacing the "off- 
diagonal" coefficients G by the truncated one G t . A similar conclusion 
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can be made for the other amplitudes in (|j). Let us estimate the error 
due to neglecting the remainder term 



G r 



sin(nfei) sm(nb 2 ) 



ly/r. 



71 



n=M+l 



rr 



(10) 



Denote 



7 := max{ 7l , 72 }, V : = ^(M+l) 2 - (N+lf 



where the maximum in the first expression is taken over the interval of 
energies we are interested in. Since N = [y/z \ , we have 

(i+M+lf-z > (l+M+l) 2 - (iV+1) 2 > {£ + V)\ 

and the relative error of the "diagonal" terms in @, i.e., the first and 
the last one, can be estimated by 



7T 



G7 



7i72 -(G*) 2 



E 

n=M+l 



-aV n 2 - 



< 



G7 



7i72-(G 



A 2 



-aD 



1-e 



-aV 



for the off-diagonal ones we have 



IA 



off 



7i72+G 2 



-aD 



-aD 



This error bounds show where deviations from the exact expression 
are most likely. An obvious requirement is that the evanescent states 
must be given opportunity to decay; as long as the moduli are of order 
of one, it is the exponential term which governs the estimates. To get 

a K digit precision, one roughly needs (M + l) 2 — (iV + l) 2 ~ (J^j . 
In general, the problem becomes therefore non-trivial when we glue 
scatterers without intermediate regions; then one has to check that the 
evanescent states do indeed decay within each Sj when we move from 
its centre to the boundaries. 

The second source of error are the denominators. The two-impurity 
system has no embedded eigenvalues (apart from the trivial ones due to 
symmetry), so the exact 7 i 7 2— G 2 is never zero in the cases of interest 
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18] . On the other hand, it has resonances unless the impurity sits at a 
node of a transverse eigenfunction; they are narrow in the case of a weak 
coupling, i.e., for a, large positive. The truncated expression may 
blow up around the resonance energies, and even if it is not the case, it 
may produce shifted resonances. This is important, because resonance 
structures such as conductivity modulations around thresholds are a 
primary object of interest in quantum- wire serial scatterers. 

To illustrate the differences due to the cut-off factorization, let us 
compute the conductivity of the quantum wire with a pair of impurities, 

o p 2 h/i] u 

g{z) = ze £ ^| W ^ ;M)|2) (n) 

11 n,m=l K n 

where t nm (z; M) is given by @ with the resolvent (|6|) which includes 
M — N evanescent states; the limit M — > oo gives the exact answer. 
Let us remark that the one-point problem contains also a sum over 
all modes in the function (||) but this can be controlled [18|. We see 



that even if the approximations with low number of evanescent modes 
i.e.M ~ 5 differ considerably only in the vicinity of the thresholds, 
there are regions like those close to a resonance (cf. the inset), where 
the convergence slows rapidly down. For instance, in order to get the 
correct position of the peak, higher evanescent modes, up to M = 15 
in this particular case, must be taken into account. 
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Figure caption: 



Figure 1: A pair of scatteres. 



Figure 2: A pair of point scatteres in the strip. 



Figure 3: Conductivity plots for factorized vs. exact scattering matri- 
ces. The positions of scatterers are ai = (0, vr/3), a 2 = (0.05,27r/3) 
and the coupling constants are chosen ot\ = a 2 = 0.2. The full line 
represents the exact result, dotted, dash-dotted and dashed lines the 
approximations with M — 1, M — 2 and M — 4, respectively. In 
the inset, however, they represent the approximations with M — 5, 
M — 10 and M = 15, respectively. 
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